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Localization of interacting fermions in the Aubry-Andre model 
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We consider interacting electrons in a one dimensional lattice with an incommensurate Aubry- 
Andre potential in the regime when the single-particle eigenstates are localized. We rigorously 
establish persistence of ground state localization in presence of weak many-body interaction, for 
almost all the chemical potentials. The proof uses a quantum many body extension of methods 
adopted for the stability of tori of nearly integrable hamiltonian systems, and relies on number- 
theoretic properties of the potential incommensurate frequency. 


Interacting fermions with a quasi-periodic potential 
describe quasicrystals or crystals with incommensurate 
charge density waves, and have a number of features in 
common with systems with random disorder; in particu¬ 
lar Anderson localization [1] can occur in the single parti¬ 
cle eigenstates. The paradigmatic model for such systems 
is the interacting Aubry-Andre model [2], which can be 
equivalently written as the Heisenberg XXZ spin chain 
with an incommensurate magnetic field. The interest 
of this model relies from one side on the fact that, de¬ 
spite its simplicity, it has a number of non trivial features 
common to more realistic systems, like a metal-insulator 
transition even at the single particle level; on the other, 
it provides an accurate description of cold atoms with 
quasi-random optical lattices generated by two or more 
incommensurate frequencies [3], [4], [5]. At the single par¬ 
ticle level, that is neglecting the many body interaction, 
its properties follow from the Harper or almost-Mathieu 
equation which have been extensively studied in the last 
thirty years; it is known in particular that the spectrum 
is a Cantor set [6] with a dense set of gaps and that its 
eigenfunctions are extended [7] or localized [8] varying the 
strength of the potential. Such remarkable properties are 
related to a deep connection between the non interacting 
Aubry-Andre model and the Kolnogorov-Arnold-Moser 
(KAM) theorem expressing the stability of invariant tori 
in quasi integrable classical hamiltonian systems. 

Much less is known when a many-body interaction is 
present, and one of the key question is which is the fate 
of Anderson localization, which is a single particle phe¬ 
nomenon. Such problem is of general interest in disor¬ 
dered quantum matter (a quasi-periodic potential can 
be considered a quasi-random disorder); in the case of 
truly random disorder analytical evidence supporting lo¬ 
calization in presence of interaction has been obtained 
[9], [10], [11], and consequences for dynamics has been 
exploited [12],[13],[14],[15],[16]; a complete proof of lo¬ 
calization in presence of interaction, even if limited to 
the ground state, is however still lacking. In the inter¬ 
acting Aubry-Andre model, analytical results have been 
established [17], [18] only in the extended regime; in the 
localized phase numerical simulations [19] and cold atoms 
experiments [5] support the stability of localization for 
attractive or repulsive interactions. 

In this letter we rigorously establish in the Aubry- 
Andre model the persistence of ground state localization 


in presence of weak many-body interaction, for almost 
all the chemical potentials. The analysis relies heavily 
on number-theoretical properties of the frequency of the 
incommensurate potential, which play a key role already 
in the non interacting case; however the structure of small 
divisors is made involved by the presence of loops caused 
by the many body interaction. Our result is obtained 
by a new technique based on a combination of multiscale 
methods developed for the classical Lindstedt series for 
KAM invariant tori in nearly integrable Hamiltonian sys¬ 
tems [20] combined with non-perturbative Renormaliza¬ 
tion Group (RG) methods [21]. 

The Hamiltonian of the interacting Aubry-Andre 
model is 

H = -e{^{a+^^aa: + ata-_-^) + fi^a+a- 

X X 

-|-C/^a+a-a-+ia;tfi (1) 

X X 

with fermionic creation or annihilation operators, 
X points on a one dimensions lattice with step 1 and 
Dirichelet boundary conditions and tpx is a potential in¬ 
commensurate with the lattice 4>x = cos(27ra;a;) and u ir¬ 
rational] it is convenient, as in the analysis of the Harper 
equation, to assume that w verifies a Diophantine condi¬ 
tion (valid on a full measure set) 

j|27rna;|j > Cn~'^ (2) 

for any integer non vanishing n ([].][ is the norm on the 
one dimensional torus). For definiteness we are choosing 
coordinates so that (f>x is even with respect to a: = 0. 

Ground state localization implies zero temperature ex¬ 
ponential decay for large values of the coordinate differ¬ 
ence in the thermodynamical correlations, for instance 
in the 2-point function < Ta“ Oy >, where T is the 
time order product, aj = , x = (a;o,a:) and 

<>= ■ In the absence of hopping or interaction 

U = e = 0 {molecular limit) there is perfect localization 
and, if < Ta“a+ > \u=e=o = 5(x,y) 

„iko{xo-yo) 

—7 - 1 - ( 3 ) 

-iko -V u(px - p- 

which is is exactly vanishing for fermions with different 
coordinates. We can now consider small U, e (compared 


5 r(x,y) = Sx,y j 


dkn 
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to e) and see if the theory is analytically close to the 
molecular limit. The interacting 2-point function can be 
written in terms of a perturbative expansion in [/, e and 
each term of the series is expressed in terms of Feyn¬ 
man graphs given by sums of products of propagators 

(3) , whose singularities naturally determine the physical 
properties. 

As the frequency w is irrational, {ujx)mod.i fills densely 
the interval (— 1 / 2 , 1 / 2 ] so that the denominator u4>x — yi 
can be arbitrarily small; this happens when {u!x)mod.i 
is close to ±(u!x)mod.i with /i = ucos27rwa;; then for 
small (wx')niod.i one has (t>x'±x — /r ~ ±(wa;')mod.i- Note 
the close analogy between the 2-point function (3) at the 
molecular limit u = e = 0, with the propagator in the free 
Fermi gas limit U = u = 0, namely {—iko+e cos k — 
the points ±a; have the same role as the Fermi momenta 
±pp such that /i = cospp, so it is natural to call them 
Fermi coordinates. This relation between molecular and 
free gas limit is a manifestation of a property of the 
Harper equation known as Aubry duality. 

We expect that the interaction causes a renormaliza¬ 
tion of the chemical potential and it is then convenient to 
fix the interacting chemical potential by choosing prop¬ 
erly its bare value; in particular the bare chemical po¬ 
tential is taken as /i = MC0s(27rwa;) -I- ev, and v is deter¬ 
mined by requiring the the dressed chemical potential is 
MCOs(27ru;a;). There are two main cases to be considered, 
corresponding to a choice of the chemical potential in a 
point of the non interacting spectrum or in a gap (which 
are a dense set). The first case is when x is irrational; in 
this case fermions living close to different Fermi points, 
that is with coordinate x' — x and x' + x with {<jjx')mod.i 
small, cannot be exactly connected by the hopping or by 
the interaction, which causes jumps of 0 , ±1 sites (but 
they can be arbitrarily close connected). Again it is con¬ 
venient to require a strong irrationality condition [7] valid 
on a full measure set, namely 

jjwa; ± 2a;x|| > C|a:|“'^ (4) 

The second case is when opposite Fermi coordinates can 
be exactly connected; this happens when 2x is an integer 
and causes the gaps in the non interacting spectrum, ex¬ 
actly like a periodic potential produces gaps in the free 
Fermi gas limit (in space of momenta instead of coordi¬ 
nate). A consequence of the above choices of x is that 
|g(A:o,a;)| < C'|x|'^, that is the denominator can be very 
small for large x; this is a sort of ultraviolet-infrared mix¬ 
ing, a phenomenon typical of incommensurate potentials. 
Our main result is the following. 

Theorem. Let us fix u = \,U = \e, and assume e, A 
small and uj verifying (2), p = (fx+s^- Then if x verifies 

(4) (chemical potentials outside the gaps of the non inter¬ 
acting spectrum) orfx is half integer (chemical potentials 
in the gaps of the non interacting spectrum), then, for a 
suitable choice ofv, the zero temperature 2-point function 
obeys to 

I < TaZa+ > I < C'e-«l"-^lF(xo - yo) (5) 


with ^ = C| logej, C a constant, and F = 

with A = (1-I-min(|x|, |j/|))“'^ if x verifies (4), while 

A = -\- and a ^ 0 if x is half integer. 

The above result establishes exponential decay in the 
2-point function, implying Anderson localization for the 
ground states of the interacting system with chemical 
potentials in the center of the gaps or in points of the 
spectrum verifying (4); in the first case there is a fast 
decay in time with rate proportional to the gap size. 

Proof. We present here the key steps of the proof (more 
technical details will be published elsewhere [17]). The 
2 -point function is obtained by the second derivative of 
the generating function 

^ j (6) 

with 

F = (7) 

e j dx(7/+V'x-l-ei +V'x-l-eiV’x) + J^ J t/” 

where if are grassmann variables, (j) is the external source, 
f dx = f dxo ei = (0,1) and P{dif) is the fermionic 
integration with propagator (3). We introduce a cut¬ 
off smooth function Xpiko,x) which is non vanishing for 
i/feg -I- {{uj{x - pS)mod.i)^) < 7> where p = ±1 and 7 > 1 
is a suitable constant (to be fixed below); therefore we 
can write the propagator as 

g{ko,x) =g^'"-'"-\ko,x) + ^ gp{ko,x) ( 8 ) 

p^± 

where gpiko,x) = and correspondingly 

tpko.x = + J2p=±i'(^P,ko,x- This simply says that 

we are rewriting the fermionic field as sum of two inde¬ 
pendent fields living close to one of the Fermi points, up 
to a regular field. We can further decompose 

0 

gp{ko,x)= gf\ko,x) (9) 

h— — oo 


with 'gp^\ko,x) similar to 'gp^\ko,x) with y replaced 
by fh with where fh{ko,x') is non vanishing in a re¬ 
gion \/kQ -\- ((wa:')modi)^ ~ 7^- After the integration 
of the generating function has the 

form 






( 10 ) 


where P{d'ip-^) has propagator g^fT^^ = Yl!k=-oo 9 p^'^ 
is given by 




ei(<h) 

pi,xop,x[+pix 


'^Pm,Xo,m,x'^+Pm 


m 


X 
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and the kernels Wm are sum of Feynman diagrams ob¬ 
tained connecting vertices e, A or ^ with propagators 
with k > h ] is given by a similar expression 

with the only difference that some of the external lines 
are associated to <j) fields. In each of the Feynman dia¬ 
grams contributing to Wm^ there is a tree of propagators 
connecting all the external lines and the vertices; the co¬ 
ordinates Xi, Xj of two external lines as such that 

* 

Xi - Xj = x'i+ PiU - x'j - PjU = '^Sa (11) 

Oc 

where the sum is over the vertices in the path of the 
tree connecting i and j and Sa = (0,1, —1) (see Fig. 1) 
is associated to the line connected to the vertex a. Of 
course Xi — Xj has an integer value. 



FIG. 1: An example of a tree; the solid lines represent prop¬ 
agators and the wiggly lines the external lines; x^i — Xw 2 = 
-f -f ^6 -f -|- ^2- 


Each kernel Wm^ is expressed by a power series in 
e,i>, X, and one has to show convergence uniformly in the 
RG step \h\; however standard power counting arguments 
says that the theory is non renormalizable, as the naive 
degree of divergence is fc — 1 if /c is the order of the graph. 
Power counting in non sufficient and one has to use a 
more accurate analysis exploiting the number-theoretic 
properties of the frequency. 

In order to do that it is convenient, given a Feynman 
graph, to consider a maximally connected subset of lines 
corresponding to propagators with scale h > hy with at 
least a scale hy, and we call it cluster v; the external lines 
have scale smaller then hy. Therefore to each Feynman 
graph is associated a hierarchy of clusters; inside each 
cluster V there are Sy maximal clusters, that is clusters 
contained only in the cluster v and not in any smaller 
one, or trivial clusters given by a single vertex. Each 
of such Sy clusters are connected by a tree of propaga¬ 
tors with scale hy; by integrating the propagators over 
the time, and using that f dxolgp^\xo, x)l < and 

that \gp^\xo,x)\ < C we get that each graph of order n 


contributing to Wm'^ 


is bounded by 


V 


( 12 ) 


where v are the clusters and hy < 0. As > 1 the 
above estimate says that the size of the graph increases 
with \h\ and diverges in the zero temperature limit. In 
getting (12), we have however not used a crucial property 
implied by the Diophantine condition: namely that if the 
denominators associated to the external lines have the 
same size, the difference of coordinates must be or 0, ±2a; 
(the second case only for x half integer) or very large; 
in this second case there is a large number of hopping 
or interactions terms by (11) and therefore a decaying 
factor associated to an high power of e. This suggests 
to distinguish in two kinds of terms, the resonant 
terms, such that the coordinate a;' measured from the 
Eermi coordinates of the external fields are all equal, and 
the remaining non resonant terms. 

We define, as usual in the theory of renormalization, a 
localization operator C acting only on the resonant terms 
and setting all the temporal coordinates of the exter¬ 
nal fields equal. We split the resonant terms in a local 
part, where C applies, and a renormalized part, where 
TZ = 1 — jC applies. In the TZ part there is at least a 
difference of fields and this 

produces a extra gain if u is a resonant cluster 

and v' is the minimal cluster containing it. Regarding 
the local terms, they are proportional to monomial of 
fields with the same coordinates x', hence anticommuta¬ 
tivity implies that terms with more than four fields are 
exactly vanishing. The possibly non vanishing quartic 
terms have the form 


sXh / dxV’1 




-■0^ , -'ll) , - (13) 

tT—,xci,X—X^ — ,Xc\,X—X V ' 


When X is irrational such terms are vanishing as there 
are fields with coordinate difference 2x which by (11) 
should be an integer. When x is half-integer local quar¬ 
tic terms exist and in such a case \h = 

Similarly the quadratic terms are of two kinds; one 
is representing the 

renormalization of the chemical potential, the other is 
d^i’X,xo,x'+pxi’+,xo,x'-px and is vanishing when 
X is irrational, while is present when x is half integer and 
in such case <Jh = 0{e'^^). 

We have now to consider the contributions from the 
non resonant terms of the effective potential In 

such terms there are least two external fields with coordi¬ 
nate a;(, ^2 with x'l ^ x' 2 ; if |h| is large the corresponding 
divisors are small but the coordinate difference is 

large. More quantitatively, if m = ^ 0^ see (11) 

and v' is the cluster containing v 

27 "-'>||(a;x;)|| + ||(a;x')||>||..(x(-x')|| = 
||w(p 2 — Pi)x + muj\\ > Co\S2x + |w||“'^ 


(14) 
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with 5 = 0 when x is irrational and 5 = 1 when is half in¬ 
teger; the last inequality follows from (2) and (4). There¬ 
fore \m\ > , that is m must be very large if the 

divisors are small and not coinciding. In a non reso¬ 
nant cluster there then are necessarily a large number of 
vertices N so that N > \m\ > by taking into 

account the hierarchy of clusters one gets 


< £2 


n 




(15) 


where are the non resonant clusters contained in v, 
that is a decay factor is associated to each non resonant 
cluster. The power counting bound (12), apparently say¬ 
ing that the theory is non renormalizable, is therefore 
replaced by 


V V V 

( 16 ) 

where ]))[* is over the resonant terms and G contains 
the product of running coupling constants. Note that 
Sv = + 5'^'^ -I- S^, where 5”^ is the number of res¬ 

onant clusters and S'® is the number of end-points with 
scale hv in the cluster v. The above bound is true not 
only for the single graph, but for the sum of all graphs 
contributing to order n and summing over all the pos¬ 
sible choices of external lines of the clusters. Note that 
[n„ 7^”'“^”] < 1. Finally one has to dis¬ 

cuss the flow of the effective running coupling constant. 
When X is irrational the only running coupling constant 
is lyfi, expressing the renormalization of the chemical po¬ 
tential and by choosing properly ly then by a hxed point 
argument lyh = When x is half integer then one 

has also a mass term and this implies 


that there is a mass scale h* of the order of the log of the 
mass. The effective coupling is = Ao(l-f 0(|h*|e) and 
as Ao = 0(e^'®“^A) then 7“^|eAh| is small. From (16) we 
have finally that the size of the n-th order contibution to 


V V 

(17) 

If 7“ >2 and e small one can then sum over {hy} and 
get an £7(7” bound. When ^ 0 this is trivial. If 
Sy^ = 0 then V is resonant and: a) when x is ratio¬ 
nal one uses the gap saying that the number of scales is 
(9(log£); b) when x is irrational, if v has two external 
lines than there is no sum as hy = hy' -I- 1 by the support 
properties of fh as the internal lines of v has the same 
ko,x' as the external ones; if has more than four lines, 
V is resonant and its local part is vanishing so one has 
an extra allowing the sum over hy. A similar 

analysis can be done for the terms with 2 external lines 
(j) contributing to the 2-point function; the exponential 
decay in the coordinate follows again from (11) as the 
only non vanishing contributions have an order greater 
than \x — 2/1/2. 

In conclusion we have rigorously established persis¬ 
tence of localization of the ground states for almost all 
the chemical potentials and in presence of attractive or 
repulsive weak interaction. The analysis relies heavily on 
the number theoretic properties of the frequency of the 
quasi-periodic potential. Our result would imply local¬ 
ization of all the eigenfunctions in a quasi-free theory, 
and we believe that an extension of the methods intro¬ 
duced here will be able to establish localization also at 
finite temperature, or extended to random disorder. 
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